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We present a theory for the Kondo spin-| effect in strongly correlated quantum dots. The theory 
is applicable at any temperature and voltage. It is based on a quadratic Keldysh effective action 
parameterized by a universal function. We provide a general analytical form for the tunneling 
density of states through this universal function for which we propose a simple microscopic model. 
We apply our theory to the highly asymmetric Anderson model with U = oo and describe its 
strong coupling hmit, weak coupling limit and crossover region within a single analytical expression. 
We compare our results with numerical renormalization group in equilibrium and with a real-time 
renormalization group out of equilibrium and show that the universal shapes of the linear and 
differential conductance obtained in our theory and in these theories are very close to each other in 
a wide range of temperatures and voltages. In particular, as in the real-time renormalization group, 
we predict that at the Kondo voltage the differential conductance is equal to 2/3 of its maximum. 

PACS numbers: 72.15.Qm, 73.63.-b, 72.10.Fk 
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The Kondo state of a quantum dot (QD) continues 
to attract significant interest from both experiment and 
theory. The zero-dimensional platform available in mod- 
ern nanotechnology provides a great advantage over bulk 
metals with magnetic impurities where the Kondo effect 
was originally observed [1] . Indeed, the possibility to en- 
rich the physics with nonequilibrium as well as to inde- 
pendently tune various parameters of QD structures al- 
lows one to access different aspects of the Kondo physics 
predicted theoretically. One of these aspects is the uni- 
versality of the zero bias anomaly or Kondo resonance 
in the differential conductance which, as has also been 
confirmed in experiments [3, l3| , turns out to be a uni- 
versal function of the temperature and voltage with the 
energy scale given by the Kondo temperature Tk. An- 
other advantage of the QD framework is the possibility to 
study the fate of the many-particle Kondo resonance in 
the presence of external fields or when the QD is coupled 
to contacts with nontrivial ground states. Recent exper- 
iments have demonstrated the universality of the Kondo 
effect also within these advanced setups, e.g., with fer- 
romagnetic contacts [J| and external magnetic fields [5|. 
Other aspects of the Kondo state, such as its use for spin 
manipulations [Gj or ferromagnetic-superconducting cor- 
relations [7[, have been addressed in modern experiments. 

The whole diversity of the Kondo physics in QDs can 
be well captured within the single-impurity Anderson 
model (SIAM) s\ . Many theoretical concepts have been 
developed to solve SIAM for QDs in the Kondo state 
in equilibrium and nonequilibrium. All these theories 
can be in general classified with respect to their appli- 
cability to two possible regimes of the Kondo state: the 
weak coupling regime, when the temperature is higher 
than the Kondo temperature, T > Tk, and the strong 
coupling regime, when T < Tk. Among the theories 
which can access only the weak coupling regime are, e.g., 
semi- analytical perturbation theories [9| well above Tk, 



or advanced analytical slave-bosonic Keldysh field inte- 
gral theories 10|, lll| which extend the applicability range 
close to Tk. In the strong coupling regime for temper- 
atures T <^ Tk there are slave-bosonic numerical mean 



field theories [12, ll3| . From high temperatures to tem- 
peratures not too much below Tk the noncrossing ap- 
proximation (NCA) [ij, ll5| is a numerical tool which, 
however, predicts a wrong scaling, Tnca t^ 2k- A less 
quantitatively reliable theory, based on the method of 
equations of motion, may be used for a qualitative de- 
scription in the whole temperature range. It has been 
applied, e.g., to QDs in an external magnetic field [16[. 
The whole temperature range can be comprehensively 
numerically described in equilibrium by the numerical 
renormalization group (NRG) method [17] which is quite 
flexible with respect to different physical setups and has 
been used, e.g., for QDs coupled to ferromagnetic con- 
tacts [18|. In nonequilibrium the recently developed 
semi- analytical real-time renormalization group (RTRG) 
method [19| is a promising theory. 

The whole spectrum of modern theoretical methods 
is in principle enough to describe many fundamental 
properties of the Kondo state in spin-i QDs in differ- 
ent regimes. However, a simple single analytical theory 
which could provide a reliable quantitative description of 
the Kondo physics at any temperature and voltage and 
which would elucidate its universality is highly desirable. 

Here we make a fundamental step towards such a the- 
ory using the Keldysh field integral formulation in the 
slave-bosonic representation. Namely, we provide the 
tunneling density of states (TDOS) in the Kondo regime, 
Eqs. dS]) and ([6]), and reduce the problem to finding a 
universal function of the temperature and voltage which 
enters this general form of the TDOS. 

We apply our theory to the highly asymmetric spin-i 
SIAM with a strong electron-electron interaction, U = 
oo (extension to finite U is straightforward), and finite 



Mo — Ed (cd is the single-particle QD energy level, no is 
the equilibrium chemical potential) at any temperature 
and voltage. The knowledge on this model up to now 
has been restricted by slave-bosonic mean field theories 
ii [H deep below Tk or, above Tk, by NCA [11 [l5| 
and the Keldysh field theory [lOJ. 

An exact description of the universal linear conduc- 
tance can be accessed within NRG. However, little at- 
tention has been paid to its prediction for the linear con- 
ductance in the highly asymmetric SIAM. But what is 
more important is that NRG fails out of equilibrium and 
the only nonequilibrium theory able to cover the whole 
voltage range, RTRG [19|, is applicable only to the s — d 
model. Thus, out of equilibrium the situation with the 
highly asymmetric SIAM with [/ = oo is even worse. 

Within the present work we eliminate this lack of 
knowledge on the Kondo regime of the highly asymmet- 
ric spin-i SIAM with C/ = oo by deriving an analytical 
expression for the TDOS which provides 1) Linear and 
differential conductance in the whole range of tempera- 
tures and voltages with the correct scaling Tk; 2) Fermi 
liquid theory (quadratic temperature and voltage depen- 
dence of the differential conductance) at low energies; 
3) Analytical ratio between the corresponding Fermi liq- 
uid coefficients. This ratio has been intensively discussed 



in the literature 20|,|21| for the symmetric spin-^ SIAM. 



However, little is known on this ratio for the highly asym- 
metric spin-i SIAM with [/ = oo; 4) Prediction identical 
to the one made in the RTRG |19j for the s — d model 
that in the deep nonequilibrium crossover region, at the 
Kondo voltage, the differential conductance is equal to 
2/3 of its maximum; 5) Excellent agreement with the 
Hamann's analytical theory 22] for the s — d model at 
high temperatures. 

Let us briefly recall the slave-bosonic approach used 
in Ref. [10|. The original QD spin-^ fermionic opera- 
tors of the Anderson Hamiltonian, d^-, d^ (u =t,i), are 
expressed in terms of new spin-i fermionic, /„, f^, and 
slave-bosonic, b, b\ operators. Since the double occu- 
pancy is forbidden, these operators satisfy the constraint, 
Q = I with Q = b^b + '^^ flfa- In practical calculations 
of observables this restriction is taken into account in the 
Keldysh field integral [23[ representation of a given QD 
observable O — T{d]y , d^) via the replacement of the QD 
Hamiltonian -ffgn with Hqd + iiQ, where /i is a positive 
real parameter with respect to which one takes the limit 
/^ — >■ c» at the end of the calculation, 

{dm = ^ hm e^^ /l?[x,x]e*^-l^"'(*)^>^"^(*)l ^^^ 

A/qM^oo J (1) 

x^r'(t),x'^(t);x^'W,x^W]- 

Here x'^^(i) and x'^i't) (and x'^'(<), x'^{t)) are the classi- 
cal and quantum components (and their conjugate part- 
ners) of the QD slave-bosonic field which is the coherent 
state field of the slave-boson annihilation operator 6, SoS 



is the Keldysh effective action which depends on /i and 
governs the dynamics of the QD coupled to the contacts, 
/? = 1/kT is the inverse temperature and A/q is the nor- 
malization constant given in Ref. 'l,(|. 

The Keldysh effective action is a nonlinear func- 
tional of the slave-bosonic fields. It represents the 
sum of the standard quadratic free slave-bosonic ac- 
tion So[x''\t),x'^{t);x''\t),x'^{t)] and the complex non- 
linear tunneling action ST[x'^^{t),x'^{t);x'^^{t)jX'^{t)] = 
~ihtTln[—iG^^^~^ — zT] , where the matrix G^") and the 
matrix T are defined in Ref. [10|: the matrix G^") de- 
scribes the isolated QD and contacts while T accounts for 
the QD-contacts tunneling coupling. The dependence of 
the tunneling action on the slave-bosonic fields x'^^i^) and 
X'i(i) comes through the matrix T. 

We now perform a time-independent shift of the clas- 
sical components of the slave-bosonic field in the matrix 

r, x'Kt) ^ x'\t) - sV2, r\t) -^ r\t) - 7V2, whue 

leaving the quantum components x'^it)^ X'^(t) unchanged. 
Note that in general 77^^ since x'^'(i) and x^^{t) are 
independent integration variables. This results in the 
appearance of nondiagonal blocks in the matrix G'"'"^ 
which now takes the form: 



G(")- 



^G!"^"\ai|(T'i') M^^\at\a'f)^ 



f(0) 



M^">{at\a't') G'c'' iat\a't') 



(o)-i. 



(2) 



.(0) 



where the inverse of the diagonal blocks, G^^, are the 
same as in Ref. [10| and the nondiagonal blocks are given 
as follows (for M\, -^ 7 is replaced with 5), 



M}^\at\at') ^ -5{t - t')-iTa, 



1 
1 



(3) 



where a is the set of quantum numbers describing the 
contacts states and Taa are the matrix elements of the 
tunneling Hamiltonian of SIAM. It is further assumed 
that a = {g,cr}, where q describes the contacts orbital 
degrees of freedom, a is the contacts s pin and Tgo-v = 
(^ct'o-t. The matrix T remains, as in Ref. [10|, off-diagonal 
but the off-diagonal blocks change (for M|, 7 is replaced 
with (5), 



Mt 



5{t-t')Ta, (r\t)-^V2 



V2h 



x'^it) 



X 



xHt) 

•\t)-l^. 



(4) 



Up to this point the theory is formally exact but obvi- 
ously cannot be solved. To make it tractable we expand 
the tunneling action St. The crucial difference with re- 
spect to Ref. [Ifll is that now we do not use the second 
order expansion in the variables x'^''''(t) (and their con- 
jugates) but instead we expand up to the second order in 
the variables x'='(t) - 5\/2, x'='(0 - 7^2 and x^)^ x\t)- 

At first sight the second order expansion of the tun- 
neling action in the variables x^^i^) ~ ^a/2; X^'^{i) — 7"\/2 
and x'^(i), x'^it) may appear inadequate with respect to 



the Kondo physics. Indeed, this expansion contains linear 
terms in the original slave-bosonic fields x'^'(t), x'^^{t) and 
X'i(i), x'^{t). These terms shift the minimum of the effec- 
tive action from the zero slave-bosonic field configuration 
to a finite one. Thus the symmetry is broken and the 
effective action of this type cannot describe the Kondo 
state having, as is well-known, no symmetry breaking. 

However, the reasoning above misses one important as- 
pect related to the Hilbert space on which the Keldysh 
effective action is defined. The point is that at fixed 
/i this Hilbert space is much wider than the QD Fock 
space having states only with zero and one electron. It is 
easy to show that when the Hilbert space is narrowed to 
the physical Fock space of the QD, i.e., when the limit 
/i -> oo is taken, the linear terms in the Keldysh effec- 
tive action do not generate any finite contribution to the 
physical observables of the QD. Thus they can be dis- 
carded from the outset and the symmetry of the Keldysh 
effective action is restored. 

In this way one obtains a Keldysh effective action 
which is quadratic with respect to the original slave- 
bosonic variables x'^'(t), x'^i't) (and their conjugates) and 
where the shifts 6 and 7 have been absorbed into the 
kernel of the action. After this transformation the the- 
ory becomes formally identical to the one in Ref. [10|. 
The only difference is that now the self-energies in the 
quadratic action are parameterized by the shifts 6 and 7. 

Therefore, one can immediately write down the general 
form for the QD TDOS in the Kondo regime, 

z..(e) = (r/27r){[6d - e + T^Kie)]' + [rj:i{e)]'r\ (5) 

where T = 27rz^c|''"P {^c is the contacts density of states). 

The crucial change with respect to Ref. [lOJ concerns 
the retarded self-energy S+(e) — SR(e) -I- iSi(e) which 
now has a nontrivial dependence on a = Sj, 



S+(e) 



■i; 



1 

2^ 



E 



ReV^ 



1 if^. W \ 

2 2TTkT 2'KkT J ' 



E^ 



^^J■x 



2TTkT 27rkT 



2TTkT / 2 



(6) 



where the sum is over the contacts {x = L,R), '4){z) the 
digamma function, /il.r = Mo T gV/2 {V is the bias volt- 
age), W is the half- width of the contacts Lorentzian den- 
sity of states and Ea = ar/2 is a complex energy. 

A simple microscopic theory for E^ = E^ + iE]^ is 
obtained using the fact that in the deep Kondo regime 
in equilibrium and at zero temperature the TDOS has 
maximum at e = /^o and this maximum is equal to the 
unitary limit 2/7rr. This leads to the two equations, 
[hi(2\E\/kTK)]siTL(t)^ {Tr/2-(t>)cos(t>, [\n{2\E\/kTK)]^ = 
7rV4 - (7r/2 - 0)^, where /cTk = 2P^ exp[-7r(^o - ed )/r] 
is the Kondo temperature, \E\ = \/(^JP"-i-T^I)^ and 
(j) = aictan{El^ / E^) . These equations show that £a = 
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FIG. 1: (Color online) The quadratic slave-bosonic Keldysh 
field integral theory (circles) provides the universal equilib- 
rium description of the Kondo spin-i effect in the asymmet- 
ric SIAM with [/ = cxD at any temperature, i.e., in the weak 
coupling regime, strong coupling regime and in the crossover 
region. The NRG results (soUd line) are provided by T. Costi 
and L. Merker (research center of Jiilich). The inset compares 
our theory with the one of Hamann [22 |. which is known to 
describe well the high temperature regime. In this plot Tk is 
defined as the temperature at which the differential conduc- 
tance maximum reaches half of its unitary limit value. One 
getsfK~1.47TK. 



Ea/kT^ is a universal ratio which guarantees the cor- 
rect scaling of the differential conductance in accordance 
with our observation above. One gets £^ ss 1.42 and 
£^ ~ 1.05. The universal ratio £a/£a provides the uni- 
versal phase, (j) ~ 2/tt. Using this model of £a we can 
calculate the differential conductance at any T and V. 

In Fig. [T] we compare the equilibrium results of our 
theory and NRG for the highly asymmetric spin-i SIAM. 
In our theory U = 00. The NRG results have been ob- 
tained for U = 250r. Increasing U further in NRG does 
not produce any significant change in the universal char- 
acteristics. In both models ^o — Cd — 8F. The pres- 
ence of the universal function in the QD TDOS leads to 
the universal temperature dependence of the differential 
conductance maximum (linear conductance) in the whole 
temperature range. As one can see from the figure, the 
universal shape of the linear conductance obtained in our 
theory is very close to the one in NRG. 

Moreover, the differential conductance as a function of 
the bias voltage turns out to be also universal. In Fig. 
[5] we show the nonequilibrium results of our theory for 
the same model as in Fig. [T]and the RTRG results [ig |. 
One sees that the nonequilibrium universal shape of the 
differential conductance for the asymmetric spin-i SIAM 
with [/ = 00 is very close to the one in the s — d model. 
In particular, in the deep nonequilibrium crossover both 
theories predict the same result for the value of the dif- 
ferential conductance at the Kondo voltage, namely, 2/3 
of its unitary limit value. This universal result has been 
used 2JI as a new experimental tool to measure Tk- 
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FIG. 2: The universal differential conductance as a function of 
the bias voltage at T = in the Keldysh field integral theory 
for the asymmetric SIAM with U = oo (circles) and in the 
RTRG theory for the s-d model (solid line). The RTRG data 
are provided by M. Pletyukhov and H. Schoeller (University 
of Aachen). The inset shows the high voltage asymptotics of 
these two models. Although the asymmetric SIAM and s — d 
model are different, their universal shapes of the differential 
conductance are very close to each other in the whole range 
of the bias voltage. In particular, both models lead to the 
identical prediction: at the Kondo voltage, eV = kTK, the 
differential conductance is equal to 2/3 of its maximum. 



At low energies we obtain (see the supplemen- 
tary material) the differential conductance analytically, 
GiT,V) = (2eV/i)[l - ct(T/Tk)2 - cvieV/kTKn The 
absence of linear terms in G{T, V) proves that the low en- 
ergy sector of our theory is the Fermi liquid. It is known 
iaim that for the symmetric SIAM cy/cr = 3/(27r2). 
However, for the highly asymmetric spin-^ SIAM with 
U = oo and finite fiQ — e^ little is known on cy/cx- For 



this asymmetric model our theory predicts 



CT 



4 21n(2|g^|) + l cv_ 

3 l^aP ' CT 



^°2^2^ 



(7) 



where c^ = [41n(2|f„|) + l]/[41n(2|£:„|) + 2] « 0.86. 

At T > Tk the linear conductance is well described 
by the analytical expression (see also the supplementary 
material) obtained by Hamann '22!]. The Hamann's the- 
ory, based on the s — d model, neglects charge fluctua- 
tions present in the SIAM. Thus it underestimates the 
conductance [25|. Our theory agrees well (inset in Fig. 



[T]) with the Hamann's result until the asymmetry and 
charge fluctuations become important at high T. 

The results above show that the nonequilibrium uni- 
versal behaviors of the highly asymmetric spin-i SIAM 
with U — oo and symmetric SIAM or s — d model are 
close to each other. The difference between them partly 
results from their fundamental physical difference {e.g., 
the Fermi liquid coefficients, like cv, are sensitive to 
the asymmetry of the SIAM; see Ref. 126[) and partly 
from the quality of our theory and the quality of RTRG. 



Within our theory the last point is equivalent to the ques- 
tion whether the theory with constant £a niay be im- 
proved through a dependence of £a on the voltage and 
temperature. In the supplementary material we show 
that a very simple dependence of f ^ on the temperature 
and voltage may, indeed, notably change the equilibrium 
and nonequilibrium results, obtained using the constant 
model. Since the function Ea = [kTKjSa effectively takes 
into account higher order terms in F, a microscopic model 
for the temperature and voltage dependence of Ea niay 
result from a proper renormalization group theory. In 
particular, taking into account the quartic terms in the 
Keldysh effective action and performing the 1- and 2-loop 
analysis may provide a renormalization of the quadratic 
term and, as a result, a function E^ = f{T/T^, eV/kT^). 

In conclusion, we have proven the existence of a uni- 
versal function S^, rigorously appearing in the formalism, 
and demonstrated how it determines the QD TDOS. A 
simple microscopic model for £„ has been proposed. This 
has allowed us to unify the strong coupling regime, weak 
coupling regime and crossover region of the Kondo state 
in equilibrium and nonequilibrium within a single an- 
alytical expression for the TDOS. To demonstrate the 
practical importance of our theory we have applied it to 
the highly asymmetric spin-i SIAM with U = oo. Our 
theory has provided the differential conductance in the 
whole range of temperatures and voltages with the cor- 
rect scaling Tk, its Fermi liquid behavior at low energies, 
an analytical expression for the ratio between the Fermi 
liquid coefficients, the prediction that at the Kondo volt- 
age the differential conductance is equal to 2/3 of its 
maximum, an excellent agreement with known theories at 
high energies. At the same time it has raised a challenge 
to develop a renormalization group method to improve 
the quality of the present theory. 
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